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Abstract 

Strong nonlinear elastic effects arise due to the Lifshitz transition associated 
with a change in the Fermi surface topology. It is shown that an electronic 
contribution to the elastic characteristics of cuprates becomes important if 
the chemical potential is near but not too near the Van Hove critical energy 
e c . In the case of a saddle point singularity in the energy spectrum an en- 
hancement factor arises for the "electronic" anharmonic constants of the 3rd 
and 4th order, that is proportional to £f/(£f — £c) to the 1st and 2nd power, 
respectively. The behaviour of the Grueneisen coefficient in the vicinity of the 
Lifshitz transition is considered. 

PACS 74.25.Ld Mechanical and acoustical properties, elasticity, and ultra- 
sonic attention, 74.25.Jb Electronic structure, 74.72.Bk Y-based cuprates, 
74.72.Dn La-based cuprates. 

Keywords: Electronic topological transition; Oxygen content; Grueneisen pa- 
rameter; High temperature superconductors. 
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High T c cuprate superconductors show unusual elastic properties in the normal as well 
as in the superconducting states. There is experimental evidence of a lattice softening in the 
normal state of La<i- x Sr : x CuO '4 (LSCO) at x — 0.14 At the same time an unexpected 
lattice stiffening has been observed in the superconducting state. For Y Ba2Cu^O x (YBCO) 
the values of the Grueneisen parameter determined by ultrasonic techniques are spread over 
an unusually broad range 0. The decrease of the Grueneisen parameter with x observed 
in [§] was recently confirmed by the oxygen concentration dependence of Raman active 
phonons || . It is of interest to consider the unusual behaviour of the Grueneisen coefficient 
as a function of the oxygen content. As we will show below, a reason for such dependence 
may be a change of electronic structure with doping. In this work we demonstrate how the 
elastic properties depend on the peculiarities of the electronic structure of the doped cuprate 
superconductors. In particular we investigate the modification of the thermodynamic and 
elastic properties due to the change of the Fermi surface (FS) topology in the nonsupercon- 
ducting metallic state. 

An essential aspect of the electronic band structure of the cuprate superconductors is 
the quasi two-dimensionality. Van Hove Q has shown that, in the 2D case, inside each band 
the electron density of states has at least one infinite logarithmic singularity (at the critical 
energy e c ). The influence of the Van Hove singularities on the properties of the high-T c 
superconductors has been widely studied ||. The effect of a Lifshitz electronic topological 
transition (LETT) in the 2D electron system of cuprates has also attracted considerable 
attention ||, [|K|, as this presents the possibility for a unified approach to studying the 
anomalies in both normal and superconducting states. The LETT occurs when the Fermi 
energy ep passes through e c , and a change in the FS topology takes place 0. The LETT 
in the 2D case for conventional metals has been investigated in J7|. The angular resolved 
photo-emission experiments || have revealed that e c is located close to Ep in the "optimally" 
doped cuprate superconductors. The possibility of the Lifshitz transition in the cuprates 
was suggested by Freeman || , and it was considered recently in the high-T c superconductors 
by means of a microscopic t — t' — J model, [H| . 



The metallic properties of cuprates are highly peculiar leading up to a deviation from 
the Fermi liquid (FL) behaviour of electrons [IT]. In particular, a dependence of magnetic 
susceptibility and Hall's constant on temperature, as well as the linear dependence of the 
resistivity on T for HTCS materials with the hole conductivity have been observed. Dif- 
ferent aspects of these anomalies have been studied in connection with the characteristic 
features of the high T c superconductors in normal phase, as it is customary to assume they 
define a behaviour of cuprates in superconducting state. There exists a pseudogap in the 
underdoped regime above T c , a rise of the localized magnetic moments and strong magnetic 
fluctuations over a wide doping range outside an area of antiferromagnetic phase. These 
effects can be considered as a manifestation of the strong electron correlations in cuprates. 
Therefore, a simple mean field treatment is not possible in this case. 

At the same time the appearance of various anomalies in the thermal and electronic 
properties of the free electron system due to the Lifshitz transition may be responsible for 
the anomalies in the interacting electron system. In this work we wish to focus our attention 
on the anomalies of the thermodynamic and elastic characteristics of cuprates which come 
from the peculiarities of the 2D electronic structure and layered anisotropy. So, we suppose 
that the FL approach is a suitable starting point as applied to normal metallic state of 
cuprates at some distance from the "optimal" doping. Below we discuss the influence of the 
electron interaction and a renormalization of the electron spectrum in cuprates. 



In the conventional metals a topological transition leads to strong nonlinear effects |12 
One can expect that elastic anomalies due to the LETT arise in the cuprate materials also. 
For simplicity, our investigation is performed for the 2D case. As well, we restrict our con- 
sideration to the case of a simple saddle point in band structure of the YBCO and LSCO 
compounds. 

To derive the behaviour of the elastic constants in the metallic normal-state of doped 
cuprates it is necessary to determine how a basic parameter of the LETT, the difference 
z = fi — e c , depends on the strain tensor Uij = 1, 2) of the 2D crystal. The renormaliza- 



tion of the chemical potential fi (at T — 0) due to deformation has been given in |T!| and it 



can be easily generalized to the 2D case. Then, we obtain for the deformed crystal 

z = z + AijUij, (1) 

— — * — * 

where zq = (sp — £ c )\u ik =o, A w - = Ay — Ay, Ay(7c) = ^ji(k) is the tensor of the deformation 
potential at tty = 0, (k is an electron wave vector), Ay denotes its average over the FS. The 
Lifshitz transition point is defined by the condition z = 0. We suppose that the term zq 
can be varied by doping so that it becomes sufficiently small. In addition, we consider an 
elastic field of deformation in the form of a sound wave. Then the transition parameter z in 



cuprates consists of a static and a dynamic part fl3fl 



The electron density of states (DOS)in the 2D system is given by the expression 

^ = 2 4/^- £ >(W' (2) 

where S is the area occupied by the electrons, 9(x) = 1 for x < 0, and else is a step- 
function. In the case of a simple saddle point, near the Van Hove critical energy e c , in the 
2D band, the electron dispersion relation may be written as 

£ ^ =£c + A (3) 

where p x ,p y are the projections of p and m\,m* 2 > are the effective masses. Substituting 
Eq. (^) in Eq. (||), it may be easily seen that a logarithmic singularity arises in v{e). We 
decompose the electron density of states into a regular part vq{s) and a singular part 8v{e): 
u(e) = fo{e) + 5u(e). The singular part results from the integration over the Fermi surface 
in the break-off region 

5v(e) = -pin(\e-e c \/e c ), (4) 

where (3 = Sn~ 2 h~ 2 (Imlm^l) 1 ^ 2 ■ For the sake of definiteness, we will consider the case when 
an equienergy surface neck appears with increasing energy, then the sign of a singular part of 
the electron density of states is positive, 5v{e) > 0. It is easy to verify that in the case of the 
neck disruption, as the energy grows, a singular part 8v{e) has the opposite sign relatively 



to Eq. (|), 8v(e) < 0. 

The free energy F can be also decomposed into a regular and a singular part: F = 
Fq + SF. 5F due to electronic transition may be represented in terms of z = p — e c , where 
p is the chemical potential, [0: 



SF = (3-\n l - 



2 Ul T n 2 T 2 



1+ 3, 2 



(5) 



here we suppose that p ^> z ^> T (temperature is in energy units). F is assumed to be 
an analytical function of z. We may neglect the temperature dependence of the elastic 
constants in the considered temperature range. 

The standard equations of the theory of elasticity are 

pd 2 u t /dt 2 = daij/dxj, a tJ = a\f + 4, (6) 

where p is a mass density, ttj is a displacement vector. The lattice term a l ^ tt holds the 
linear contribution Xiji m Uij, where \ji m is a tensor of elastic moduli, and the conventional 
anharmonic terms of expansion. The electronic term a[ ~ d5F/ duij contains a logarithmic 
singularity, 



a'ij = (3 zAij In — + ^zA i;j . (7) 



The constant (3q = (3 /V is of the order ~ vo{ep) /tt 2 Vq, where u (e_p) is the regular part of 
density on the FS, Vq is the volume of the unit cell of the crystal. 

There are two different limiting cases from Eq. ([I]): (i) \zq\ = \ef — £ c \u i:j =o is comparable 
with |Ajjitjj| in magnitude, or (ii) far exceeds it. We consider the latter case, when the sound 
wave propagates along the metallic phase of the cuprates, and the values Uij are such that 

\z \ > \AijUij\- (8) 

Then, using the expansion of SF, Eq. (|5|), with respect to the small parameter | A^-w^- 1 / 1 ^ 1 ^ 
1, (here the singular terms at \zo\/e c ^ 1 are retained) 

SF - -(AijUij) 2 In M ± ^ K ^ m Y _ HA lm ui m ) 4 { ) 

2 V 3 3! e c 6 \zq\ 24 z 2 ' v ; 



the electronic stress tensor, Eq. (|7|), may be expanded up to the second order terms as: 

o\ • = PoAijAirnUim In — ± ^ ^-(A lm Ui m ) 2 . (10) 



Here the upper sign corresponds to the region zq > where the Fermi surface neck exists, and 
a lower sign refers to region z$ < where the neck has been disrupted. The quadratic term 
m Eq. (0) has the same form as the lattice contribution due to the cubic anharmonicity. 
Using Eq. (||) and the estimate \Aij\ ~ A5ij, where A > is of the order of the deformation 
potential constant, the "electronic" anharmonic constants of 3rd and 4th order, 1^, 1^, can 
be written as 

r 3 ^±^A 2 M, r 4 --^AV, (ii) 

6 24 y J 

where \r)\ = A/ \z \ ~ epj \ep — e c \ 3> 1; A > e F as we shall discuss below. The case 
considered here could be formally interpreted as an increase of the anharmonic constants 
T 3 , T 4 with an enhancement factor rj, to the corresponding power. Thus the electronic 
contribution to the anharmonicity increases as the LETT point is approached. 

In conventional metals A is usually of the same order as Bp [|12|| . Since a value of Bp in 



cuprates is well below, at first sight it would seem that the influence of electrons on elastic 
properties of these materials is modest. But in the metallic state of cuprates, owing to a 
relatively low density of carriers and a poor screening (compared to the standard metals), 
the Coulomb energy of interaction is large, and can cause the value of A to increase. The 
quasi-two dimensionality of the electronic structure in cuprate superconductors also provides 
an explanation for A increasing. For the quantitative estimates of A we used a familiar 



procedure from the method of deformation potential ||14j| . Hence we take into account that 
the electrons of conductivity in cuprates are basically localized in the C11O2 planes. Using 
an effective density of electronic states per atom of Cu, vq(ep) = 4.8 (eV • at.Cu) -1 in 
the case of YBCO and v (e F ) ^ 4.14 (eV • at.Cu)- 1 for LSCO, |TJ, we obtain A =1.25 
eV and A = 0.82 eV, respectively, which are larger than the accepted values of Ep in these 
systems. 



The relation between the energy distance zq and the doping parameter x is of great 



importance ]nj. For simplicity, we treat the dependence of Zq on the doping x as an 
empirical fact. The dopant dependence of zq has been estimated from experimental data for 



the chemical potential shift in doped La2- x Sr x CuO^ ]T7|], and Y Ba2Cu^O X} [|T^]. From 
the above-mentioned estimates we have obtained the following values of the "electronic" 
anharmonic constants in LSCO for \z Q \ ~ 100K: |r 3 | = 1.33GPa, T 4 = -30GPa; in YBCO 
for \zq\ ~ 250K: | T3 1 = 14GPa, = — 189GPa. The dependence of "electronic" anharmonic 
elastic constants of 3rd and 4th order for compounds LSCO and YBCO on doping are 
shown in Fig.l (b), (d). A dependence of T^, on doping, Eq. (|i~T|) , and the asymmetry of 
T 3 relatively to the point Zq = offer a reason for a distinction between a metallic phase of 
cuprates and the conventional metals in the kinetic characteristics and transport properties. 



Also, the anharmonic dynamics is highly sensitive to the ion mass fll5| , so that one can 
expect a modification of isotope shift with doping. 

The Grueneisen constant is defined as the ratio of the full coefficient of thermal expansion 
ay to the full specific heat of the crystal at constant volume Cy 

VayB T V ( d 2 F \ d 2 F 

7 = ^-, or 7 = _^_j/_ (12) 

where B? is the isothermal bulk modulus and F is the free energy. The spectrum of oscilla- 
tions in cuprates as anisotropic layered crystal structure is characterized by more that one 
Debye temperature. In particular, a temperature range, where the interaction between the 
layers is not essential, is |19| , 

(£ ± /£|,) 1/2 6<T<0, £±<% (13) 

here @ = nTis\\/a is a characteristic temperature, corresponding to a velocity sn of propa- 
gation of sound wave in the afr-plane of layer C11O2, s\\ = {E\\/ p) 1 ^ 2 ; a is a cell parameter, 
E\\ and E± are the Young's moduli in the plane ab (X xxxx ) and along the c-axis (X zzzz ) 
accordingly. Let us estimate the temperature range defined by Eq. (|T3"D. In the case of 
YBa 2 Cu 3 7 , using Eh = 2.30 x 10 n Pa, E ± = 1.60 x 10 n Pa E0|, and a mass density 



p = 6.4g/cm 3 , we obtain 6 = 356K and (E ± /E ]l ) 1 / 2 Q = 296K. For LSCO respectively, 
using E\\ = 2.48 x 10 n Pa, E ± = 2.05 x 10 n Pa ffift, p = 6.88g/cm 3 , we find = 361K 
and (E±/E\\) 1 / 2 = 328K. Notice that the temperature boundaries in Eq. (|L3) in turn 
themselves depend on the doping. 

A singular part of the lattice Grueneisen parameter <57 p?t can be obtained by estimating 
the corrections to the squared phonon frequency averaged over the Brillouin zone, resulting 
from the topological transition. Our estimates have shown that the latter are relatively 
small. Much more sensitivity to the doping in the region of the LETT was found for the 
electronic Grueneisen coefficient, r y e . To calculate 7 e , we shall present, as previously, the 
free energy as a sum F = E + SF, then 7 e = 7q + ^7 e - 

The regular parts of the electronic specific heat Cy and the thermal coefficient of pressure 



arc 



According to Eq. ([5]), the singular part of the electronic specific heat follows as 

^bF ^ ir 2 (3 \ v \e c 
SC V = -T— = -—Tin— (15) 

The singular part of the thermal coefficient of pressure is 

lJdP\ 1 d 2 6F „ vr 2 /? 

T 5 \&T) v = -TWdV = -— ri - (16) 

Using Eqs. (0), (@)- (0) and detaching a standard electronic Grueneisen coefficient 7q = 
~jjy = d\niy /d\nV , we obtain the electronic contribution to the Grueneisen parameter 
near the point of the LETT as 

7 e = l e /c(v) + ^7 e , 57 e = -klA 2 c(^), (17) 

where c(r/) = (1 + ^ In ^)(1 + Cf/Cy). The dependence of the singular electronic cor- 
rection 5Y due to the proximity to the LETT on doping accordingly to Eq. (|T7D is plotted 
in Fig.l (a), (c). Under certain conditions its variations may lead to negative values of the 



Grueneisen coefficient. In particular, important point is the relation between the phonon and 
electronic specific heats C v h /C v which enters in Eq. (|TT|). At relatively high temperatures 
(T > T c ) the phonon contribution to the specific heat C v is, as a rule, greater than the elec- 
tronic one, C v . But in the temperature range Eq. (|T3|), where C v h ~ T/Q due to the layered 
anisotropy C v h is comparable with C v . For the metallic normal state of LSCO another 
low temperature interval where C v 3> Cy 1 is available: T <C (e±/E\\ \ ^ 6(6 / 'ep) 1 ^ 2 '■ This 
case is the only one plotted in Fig.l (a). 

Let us take up the assumptions underlying the results obtained above. We assumed 
Bp \z\ ^> T. The restriction ep ^> \z\ is obvious from setting up a problem, since the 
ratio |z|/e_F is a dimensionless parameter that defines the distance to the transition point. 
From Eq. (|j) it follows that \uij\ <C |z |/|Ay| ~ ko|/A <C 1. Thus, the present treatment 
is appropriate to ultrasound waves (with finite amplitudes) as well as to long-wavelength 
phonons. In this connection of special interest is a problem of interplay between the "elec- 
tronic" anharmonicity in plane ab due to the proximity to the Lifshitz transition and a 
"tilting" mode in LSCO, or, in case of YBCO, a pyramidal apex oxygen motion along the 
c-axis. 

For temperatures in considered T <^ Q D , z & boundary of superconducting state gives a 
restriction from below, T > T c . At T 7^ the singularities caused by the LETT are smeared 
out. The width of the temperature smearing is Az ~ T. So, the obtained results are valid 
at temperatures such that the smearing is small: T/|^o| ~ ^z/\zq\ 1. The inequality 
T/\zq\ 1 coincides with a condition for a degeneracy of the electron gas in the normal 
metallic state of cuprates. The Fermi-liquid approach in this case may be considered as 
a close approximation if A (6d/T) 2 f22f , where A is an electron-phonon coupling 
constant. The simple estimates of A for YBCO and LSCO show that above mentioned re- 
quirement is fulfilled for these compounds at A ~ 1.3 -j- 4.2. While a behaviour of A depends 
on the details of anharmonic interactions ||23|| . Taking into account the "electronic" contri- 
bution to the anharmonic constants, Eq. (|TJ), one can find A as a function of temperature 



24j[ and doping, to refine the regions of the Fermi-liquid theory feasibility. 



Up to now we considered a free electron system in the presence of the simple saddle point 
in the electron spectrum. Now we discuss how the fact that the electrons in the cuprates are 
a strongly correlated system will influence the results. Recently in [f25j it was shown that 
in 2D interacting system the density of states deviates from the bare DOS. Particularly, for 
the doping range around the critical point corresponding to the LETT, it is characterised 
by a shifted logarithmic singularity and two jumps on the right in energy of the logarithmic 
divergence point. Ignoring these jumps [], one can conclude that the obtained results are 
qualitatively not changed for a renormalized spectrum in the presence of interaction Q. The 
quantitative changes are appreciably defined by the magnitude of the deformation potential. 
At the same time the renormalization of the deformation potential in the system of strong 
correlated electrons remains to be solved and demands further investigation. In this casen, 
it is of special interest to extend the deformation potential theorem for the electron phonon 
coupling due to Bardin to the electron-spin fluctuation coupling as given by Schrieffer |2(| . 

The comparison of the anharmonic constants of 3rd and 4th order, r 3 ,r 4 , Eq.(ll), 
with the experimental values is not possible because of lack of data. However there are 
indications that the results obtained with the use of T 3 , T 4 are consistent with other experi- 
mental data. Fig.l (c) shows that a change of the singular correction, £7% of the electronic 
Grueneisen coefficient in Y Ba2Cu^O x , as function of oxygen content x (solid curve) is in 
qualitative agreement with a variation of a mean acoustic Grueneisen parameter with x 
taken from experimental data 0. A closer comparison would require further experimental 
data with monocrystal samples in the proximity to the "optimal" doping. In the case of 
La,2- x Sr x CuOi ) at x=0.14 the experimentally observed anomaly (negative value) of the in- 
plane ab thermal expansion coefficient an |27j gives one more evidence of agreement between 



our results and experiment. When ocu is expressed in terms of the Grueneisen parameter 



1 What means that we have restricted ourselves to the region of the overdoped regime. 

2 A role of the jumps in the underdoped regime for elastic properties we consider elsewhere. 
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using Eqs. flT2|), (p7|) , we obtain ay < in the required temperature range. 

In summary, we have presented the dependence of the elastic properties on the singular- 
ities of the 2D electronic energy spectrum in the metallic normal state of the doped cuprate 
superconductors, as applied to LSCO and YBCO compounds. It has been shown that the 
proximity to the point of a topological modification of the Fermi surface gives rise to the un- 
usual "electronic" anharmonicity, which essentially depends on the doping and may exceed 
the usual lattice anharmonicity. The correction to the electronic Grueneisen coefficient due 
to this anomalous anharmonicity was found to be negative, what is in line with the exper- 
imental data for YBCO and LSCO. In superconducting state, "electronic" anharmonicity 
can be of importance in the context of an enhancement of the electron-phonon coupling to 
explain high T c in cuprates. 
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Grant no.436-Rus-113/202/0/R,S. V.R. gratefully acknowledges the hospitality of the Uni- 
versity of Rostock. 
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FIGURES 
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FIG. 1. Dependence of singular correction 5^f e on doping x for LSCO (a) and for YBCO (c) 
near the Lifshitz transition point. In (b), (d) the similar dependence of elastic constants I^, T± is 
shown. Dashed lines separate the regions with zq > (left), zq < (right). In (c) the open circles 
are the experimental points for a change of mean acoustic Grueneisen parameter with x 0. A 
dotted line is intended as a guide to the eye. 
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